Mesoscopic Stoner instability in metallic nanoparticles revealed by shot noise 
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We study sequential tunneling through a metallic nanoparticle close to the Stoner instability 
coupled to parallely magnetized electrodes. Increasing the bias voltage successively opens transport 
channels associated with excitations of the nanoparticle's total spin. For the current this leads 
just to a steplike increase. The Fano factor, in contrast, shows oscillations between large super- 
Poissonian and sub-Poissonian values as a function of bias voltage. We explain the enhanced Fano 
factor in terms of generalized random-telegraph noise and propose the shot noise as a convenient 
tool to probe the mesoscopic Stoner instability. 
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Introduction. - The onset of ferromagnetic long-range 
order in bulk metals was proposed by Stoner [1] as a 
consequence of the competition between kinetic and ex- 
change energy. The Stoner instability leads to a state 
with a macroscopically large total spin S. This may be 
different in mesoscopic systems such as metallic nanopar- 
ticles or large quantum dots, which can be modeled 
by a universal Hamiltonian characterizing their discrete 
single-particle level spacing. Coulomb interaction, and 
exchange interaction of the electron spins. The latter 
gives rise to the so called mesoscopic Stoner instabil- 
ity. The nanoparticle becomes partially polarized and 
the value of the total spin S sensitively depends on the 
ratio of exchange interaction and level spacing [2, 3]. Re- 
cent theoretical studies of the mesoscopic Stoner regime 
addressed the statistics of linear conductance peaks [4], 
the interplay with Coulomb charging [5, 6] and the Kondo 
effect [7, 8], the influence of spin-orbit interaction [9, 10], 
the competition with superconducting pairing [11, 12], its 
sensitivity to the degree of chaos [13] and its implications 
on the thermopower [14]. 

An experimental confirmation of the mesoscopic 
Stoner instability is difficult for at least two reasons. 
First, magnetic moments of individual nanoparticles or 
quantum dots are hard to measure. Second, the exchange 
coupling of a given nanoparticle is not tunable, which 
would be desirable for studying the successive transitions 
from zero to full polarization. So far, only the behavior 
of the Coulomb peaks as a function of magnetic field 
served as a tool to probe the spin properties of quantum 
dots [15-18]. 

Transitions between different values S of the total spin 
can also be achieved by tunnel coupling the nanoparticle 
to source and drain electrodes with a finite bias voltage 
V . Increasing the bias voltage enables transitions be- 
tween different values S of the total spin. In this Letter, 
we suggest to identify the mesoscopic Stoner instability 
by using ferromagnetic electrodes with parallel magneti- 
zation directions and measuring the shot noise. We find 



that for an even number of accessible values for the to- 
tal spin 5, the Fano factor is strongly enhanced while 
for an odd number it stays sub-Poissonian. This results 
in oscillations of the Fano factor as a function of bias 
voltage that are robust against asymmetries in the tun- 
nel couplings. We propose these oscillations as a striking 
evidence of the mesoscopic Stoner instability. 

Model and technique.- The nanoparticle is repre- 
sented as a multi-level quantum dot described by the 
universal Hamiltonian [2, 3], 



Hdot = eacLw + Ec{N - Naf 



(1) 



The first term describes the spin-degenerate single- 
particle levels a of the quantum dot. The second term 
models Coulomb charging. Here, Ec is the scale of the 
charging energy, N = J2aa '^aa'^acr is the number of elec- 
trons on the dot, and Nq is the equilibrium charge of 
the dot that can be tuned by a gate voltage. Finally, 
the third term describes a Heisenberg-type spin interac- 
tion for the total dot spin S = J2aaa' c\,aCraa'Cacr' with 
ferromagnetic exchange coupling J. 

In the following, we assume a constant single-particle 
level spacing A and the following hierarchy of energy 
scales: Ec,A,J > > A - J > fceT. The bias 
voltage V defines the energy scale for the available exci- 
tations. The inequality Ec 3> eV implies that only two 
charge states A'o and A'o + 1 participate in transport. 
Due to the condition A ^ eV we can neglect all states 
with particle-hole excitations. Hence, the dot eigenstates 
|A^, S, Sz) are fully characterized by the total number N 
of electrons, the total spin S and its z-component Sz. 
The corresponding eigenenergy. 
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is independent of Sz. For a given number of electrons 
iV, the ground state spin S is the integer value (for even 
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N) or half-integer (for odd A'') value that is closest to 
J/[2(A — J)]. This value increases with increasing J and 
diverges for J ^ A. The latter marks the (macroscopic) 
Stoner instability. The excitation energies for changing 
the total spin S by one while keeping the charge number 
A'' fixed form a ladder with constant spacing 2(A — J), 
i.e., the condition eV > A — J ensures that more and 
more spin states are involved with increasing bias voltage. 
Finally, the inequality A — J ^ k^T guarantees that 
thermal smearing is small enough to resolve individual 
states with different total spin S. 

To achieve transport, the nanoparticle is tunnel cou- 
pled to two ferromagnetic electrodes with parallel mag- 
netizations. The Hamiltonian of the total system is given 
by H = Hdot + J2r + ^^tun- The ferromagnets are de- 
scribed in terms of noninteracting electrons at chemical 
potential as Hr = Sko- ^k<Tarkcr"'-k<T, where al^^ cre- 
ates an electron with momentum k and spin a in lead 
r = L,R. We assume the leads to have a constant, spin- 
dcpcndcnt density of states pra which is related to the 
polarization of the leads via = {Pr+—Pr-)/{Pr++Pr-)- 
In the following, we assume pl = Pr = P- The cou- 
pling between the dot and the leads is described by 



-fftun = J2rkaa'tral^^^Cacr + H.c. The tuuncl matrix el- 
ements tr are related to the tunnel-coupling strength 
via Tra- = 27r I Pro- ■ Furthermore, we define F^ — 

E.r../2. 

In order to evaluate the current and the current noise 
we use a real-time diagrammatic approach [19, 20] which 
takes into account the interactions on the dot exactly and 
performs a perturbative expansion in the dot-lead cou- 
pling. The idea of this approach is to integrate out the 
noninteracting lead electrons and to describe the remain- 
ing quantum-dot system in terms of its reduced density 
matrix. It is sufficient to consider only diagonal matrix 
elements Pn,s,S^ which describe the probability to find 
the system in the eigenstate \N,S,Sz) [21]. In the sta- 
tionary state, they are determined by the master equa- 
tion 

0= E W^''fN--''''pN,s,s., (3) 

where Wy^ are the transition rates from \x) = \N,S, Sz) 
to Ix') = l-^'j 5", S'z) evaluated to first order in the tunnel 
couplings F^. The rates for an electron entering the dot 
are given by 
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where (5', 5'^; i, ±i|S", S'^) denotes the Clebsch-Gordan 
coefficient for adding the spin of the incoming spin 
up/down electron to the initial spin S, Sz to get the final 
spin S', S'z, and fr{x) = l/(e(^-'''-)/'=BT + i) ig the Fermi 
function. The electrochemical potentials of the leads are 
chosen symmetrically, /iR — — /iL = eV/2, such that elec- 
trons preferably travel from the left to the right lead. 

The rate W^'.'^^i for the reverse process is the same 
S' S' 'S s 

as WJ^^-^^'J^' ^ but with fr being replaced by 1 — fr- 

Prom the form of Eq. (4), that applies for parallel 
magnetized electrodes, it follows that the probabilities 
Pn,s,S^ are independent of Sz, i.e., the spin state of the 
nanoparticle remains rotationally invariant even under 
voltage bias. A different relative orientation of the leads' 
magnetizations, e.g., antiparallel, would break this sym- 
metry and allow for accumulation of a finite dipole mo- 
ment. 

The current through the system is given by / = (/) = 

/ s' S' ' s s 

Y.N,N',s,S',s^,S',^ n']n' ' 'PN,s,S,, where the current 
rates are obtained from the transition rates by multi- 



plying with the number of electrons transferred between 
the dot and leads. We use standard techniques to cal- 
culate the current noise S = J dt{i{t)i{0) + /(0)/(i) - 
2(/)^) in Coulomb-blockade systems for weak tunnel cou- 
pling [22-24] which may also be formulated within a dia- 
grammatic language [25]. From this we obtain the Fano 
factor defined as F = S/{2eI). 

Results.- For large charging energy Ec, only two 
charge states, A^ and A' -I- 1, are involved in transport. 
The number M of possible values for the total spin S 
increases with increasing bias voltage. An example for 
the available states and how they are connected by tun- 
neling is sketched in Fig. 1. The interesting regime for 
detecting the mesoscopic Stoner instability is achieved 
when Ng is tuned (via a gate voltage) in such a way 
that the lowest excitation energy for a charge transfer, 
e.g., -EAr+i,s+i/2 — En^s, is smaller than the energy for 
spin excitations without charge transfer, En^s±i — -E-Af.Sj 
where A^ and S are the ground state charge and total 
spin in the absence of bias voltage. In this case, M in- 
creases successively one by one with increasing bias volt- 
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FIG. 1. (Color online) Nanoparticle states and possible tran- 
sitions. Solid (dotted) lines indicate tunneling of majority 
(minority) spin electrons. 




FIG. 2. (Color online) Current as a function of bias voltage. 
Arrows mark onset of transport through excited states. No = 
n+ 1 + ^^;^ -0.0001, n€'N,Ec = lOA, J = 0.97A, Fl = 
Fr, p = 0.3, T — O.OOIA. Higher steps become increasingly 
difficult to resolve. 



age. This is accompanied with a step in the current, see 
Fig. 2, which is smeared by temperature. Such a series of 
steps occurs whenever excited states, e.g., of vibrational 
or magnonic nature [26, 27], contribute to transport and 
is no unique signature of the mesoscopic Stoner regime. 

Much more interesting behavior is seen in the Fano 
factor, see Fig. 3. The latter oscillates between large 
super-Poissonian (for even M) and sub-Poissonian (for 
odd M) values. The (zero-temperature) Fano factor of 
the first plateau {M = 2), can be written in a compact 
analytic form. For the case that only the states | A^, S, S^) 
and |Af + 1, S" + ^, Sz) contribute to transport, we find 
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where 7l = (25 -I- 2)Fl and 7r = (25 + 1)Fr. 

In the limit of unpolarized electrodes, p = 0, the Fano 
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FIG. 3. (Color online) Fano factor as a function of bias volt- 
age for different asymmetries of the tunnel couplings. Other 
parameters as in Fig. 2. 



factor F = (7l +7r)/(7l +7r,)^ remains sub-Poissonian. 
For the limits 5 = and 5 — > oo, we recover the well- 
known results for transport through a single-level quan- 
tum dot, F = (4F2 -|-F|)/(2Fl + Fr)2 [28] and the resuh 
for spinless electrons, F = (F^ -|-F2,)/(Fl +Fr)2 [29], re- 
spectively. 

To achieve a super-Poissonian Fano factor, a finite po- 
larization, p 7^ 0, is required. For 5 = we recover 
the result known for transport through a single-level 
quantum dot, F = [4(1 + p^)Tl + {I - p^)TW/[{l - 
P^)(2Fl + Pr)^] [30], which becomes super-Poissonian 
for p > p* = ^Fr/ (2Fl + Fr) and even diverges for 
p — > 1. The super-Poissonian noise (for 5 = 0) is caused 
by bunching of electrons [31, 32]. Majority-spin electrons 
are transferred through the dot until finally a minority- 
spin electron enters. It blocks further transport as it has 
a small probability to leave to the drain. Therefore, the 
current gets chopped into bunches of majority-spin elec- 
trons. 

For probing the mesocopic Stoner instability, we con- 
sider the case of finite p and large 5. The leading contri- 
bution 
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F= -^-^52 + 0(5) 
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to the Fano factor (of the first plateau) scales quadrati- 
cally with 5, diverges for p — > 1, and is independent of 
the tunnel couplings. This result can be understood as 
a generalization of random-telegraph noise. For this, we 
consider the different projections Sz for given N and 5 as 
the different states among which the system can switch. 
As remarked above, the probability distribution pn,s,S^ 
is 52-independent. Current is established by transitions 
that change the state to\N + 1,S + 1/2, Sz ± 1/2). For 
p = 0, the two possibilities to increase or decrease the 
z-component of the spin, weighted by the corresponding 
Clebsch-Gordan coefficients, add up such that for each 
of the states not only the probability but also the con- 
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tribution to the current is independent of the initial Sz- 
This is different for finite polarization, p ^ 0. Then, 
transitions involving majority-spin electrons (indicated 
by solid lines in Fig. 1) are more likely than those involv- 
ing minority spins (dotted lines). As a consequence, the 
contribution to the current depends on the initial value 
of Sz- In the considered example, it monotonically in- 
creases with Sz- For a qualitative estimate of the Fano 
factor, we assume symmetric couplings, Fl — Fr , neglect 
the Clebsch-Gordan coefficients, and perform a mapping 
onto an effective two-state system. The state represent- 
ing the positive/negative values of Sz carries a current 
I± = e(l ±p/2)Fl. The average time to change Sz by 
1/2 is To = 1/[(1 - p^)Fl]. The time r to switch from 
positive to negative Sz by random walk scales with the 
square of the distance, i.e., r = {2S)^to. For frequencies 
smaller than this switching time, the noise of a random- 
telegraph signal is given by (/+ — /_ )^r/2 [33] from which 
we estimate the Fano factor as 5*^/(1 — p^) in qualita- 
tive agreement with the result in Eq. (6). 

With increasing bias voltage, a third value of the total 
spin, e.g., 5—1/2 in addition to S and S + 1/2, becomes 
available (M = 3) and the Fano factor drops to some sub- 
Poissonian value, see Fig. 3. From the alternating way 
of how majority and minority spins are involved for the 
transitions it follows that the contribution to the current 
monotonically increases with Sz for total spin 5'+ 1/2 but 
decreases for S — 1/2, while it is roughly independent of 
Sz for S. As a consequence, the system can rather quickly 
change between the high- and low-current states and the 
random-telegraph noise is reduced. 

For larger values of M, an S'^-dependence of the cur- 
rent is restricted to the smallest (S'min) and largest 
{Smax — Sn^in + [M — l)/2) valuc of the total spin. A 
super-Poissonian Fano factor appears for even M, since 
the current either increases or decreases with Sz for both 
S'min and Smax- In coutrast, the Fano factor remains sub- 
Poissonian for odd M, since the current increases with Sz 
for S'min but decreases for Smax, or vice versa. 

We finally address the effects of asymmetric tunnel 
couplings on our results. To this end, we show in Fig. 3 
the Fano factor as a function of the applied bias voltage 
for different asymmetries of the tunnel couplings. Inter- 
estingly, the oscillations of the Fano factor persist even 
for large asymmetries. The asymmetry only has the ef- 
fect of increasing the Fano factor slightly when many 
spin states contribute to transport. Thus, the predicted 
Fano factor oscillations turn out to be very robust which 
should facilitate their experimental observation. 

Conclusions.- In conclusion, we discussed spin- 
polarized transport through a metallic nanoparticle close 
to the mesoscopic Stoner instability. When only two spin 
states contribute to transport, we found that random- 
telegraph noise leads to a super-Poissonian Fano factor 
that scales with the square of the total spin. We, further- 
more, found that the Fano factor oscillates as a function 



of bias voltage showing super (sub-)Poissonian behavior 
when an even (odd) number of spin states contributes to 
transport. These oscillations turn out to be robust with 
respect to asymmetries of the tunnel coupling, thereby of- 
fering a good tool to experimentally investigate the meso- 
scopic Stoner instability. 
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